Orbital Fluctuation Mediated Superconductivity in Iron Pnictides: 
Analysis of Five Orbital Hubbard-Holstein Model 
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In iron pnictides, we find that the moderate electron-phonon interaction due to the Fe-ion os- 
cillation can induce the critical d-orbital fluctuations, without being prohibited by the Coulomb 
interaction. These fluctuations give rise to the strong pairing interaction for the s-wave supercon- 
ducting (SC) state without sign reversal (s++-wave state), which is consistent with experimentally 
observed robustness of superconductivity against impurities. When the magnetic fluctuations due 
to Coulomb interaction are also strong, the SC state shows a smooth crossover from the s-wave 
state with sign reversal (s±-wave state) to the s++-wave state as impurity concentration increases. 

PACS numbers: 74.70.Xa, 74.20.-z, 74.20.Rp 



The mechanism of high-T c superconductivity in iron 
pnictides has been an important open problem. By con- 
sidering the Coulomb interaction at Fe-ions, antiferro- 
magnetic (AFM) fluctuation mediated fully-gapped sign- 
reversing s-wave state (s±-wave state) is expected theo- 
retically [1, 2]. Regardless of the beauty of the mecha- 
nism, there are several serious discrepancies for the s±- 
wave state. For example, although s±-wave state is ex- 
pected to be very fragile against impurities due to the 
interband scattering [3], the superconducting (SC) state 
is remarkably robust against impurities [4] and a-particle 
irradiation [5]. Moreover, clear "resonance-like" peak 
structure observed by neutron scattering measurements 
[6] is reproduced by considering the strong correlation ef- 
fect via quasiparticle damping, without the necessity of 
sign reversal in the SC gap [7] . These facts indicate that 
a conventional s-wave state without sign reversal (s++- 
wave state) is also a possible candidate for iron pnictides. 

Then, a natural question is whether the electron- 
phonon (e-ph) interaction is important or not. Although 
first principle study predicts small e-ph coupling constant 
A ~ 0.21 [8], several experiments indicate the significance 
of e-ph interaction. For example, the structural transi- 
tion temperature T§ is higher than the Neel temperature 
in underdoped compounds, although the structural dis- 
tortion is small. Also, prominent softening of shear mod- 
ulus is observed towards Ts or T c in Bal22 [9]. Raman 
spectroscopy [10] also indicates larger e-ph interaction. 

Interestingly, there are several "high-T c " compounds 
with nodal SC gap structure, like BaFe 2 (Asi_ a; P, E )2 [11] 
and some 122 systems [12]. Although nodal s±-wave 
state can appear in the spin-fluctuation scenario due to 
the competition between the dominant Q = (it, 0) and 
subdominant fluctuations [1, 13], the T c is predicted to 
be very low. Thus, it is a crucial challenge to explain the 
rich variety of the gap structure in high-Tc compounds. 

In this letter, we introduce the five-orbital Hubbard- 
Holstein (HH) model for iron pnictides, considering the 
e-ph interaction by Fe-ion vibrations. We reveal that a 



relatively small e-ph interaction (A < 0.3) induces the 
large orbital fluctuations, which can realize the high-T c 
s ++ -wave SC state. Moreover, the orbital fluctuations 
are accelerated by Coulomb interaction. In the presence 
of impurities, the s++-wave state dominates the s±-wave 
state for wide range of parameters. 

First, we derive the e-ph iteration term, consider- 
ing only Einstein-type Fe-ion oscillations for simplic- 
ity. Here, we describe the <i-orbitals in the XYZ- 
coordinate [1], which is rotated by 7r/4 from the xyz- 
coordinate given by the Fe-site square lattice: We write 
Z 2 , XZ, YZ, X 2 -Y 2 , and XY orbitals as 1, 2, 3, 4, 
and 5, respectively [1]. We calculate the e-ph matrix 
elements due to the Coulomb potential, by following 
Ref. [14]. The potential for a d-electron at r (with 
the origin at the center of Fe-ion) due to the surround- 
ing As 3 ~-ion tetrahedron is [7 ± (r;tt) = 3e 2 J2t=i \ r + 
u — Rfl^ 1 , where u is the displacement vector of the 
Fe-ion, and Rf is the location of surrounding As-ions; 
V3R+/R F c-As = (±V2,0,1) and (0,±\/2,-l) for FeW, 
and VSR7/Rfc-a s = (±%/2,0,-l) and (0,±%/2,l) for 
Fe^ 2 ) in the unit cell with two Fe-sites. Note that ux,y 
and uz belong to E g and B\ g phonons [10]. The u linear 
term of , which gives the e-ph interaction, is obtained 
as V ± (r; u) — ±A[2XZ-u x — 2YZ ■ u Y + (X 2 — Y 2 )u z ] + 
0(r 4 ), where A = 30e 2 /v^f4 c _ As . Then, its nonzero 
matrix elements are given as 

(2|V|4) = ±2a 2 Au x /7, (3|F|4) = ±2a 2 Au Y /7, 
(2\V\2) = ±2a 2 Au z /7, (3\V\3) = t2cl 2 Au z /7, (1) 

where a is the radius of <i-orbital. Here, we consider 
only for orbitals i, j = 2 ~ 4 that compose the 

Fermi surfaces (FSs) in Fig. 1 (a) [1]. The obtained e-ph 

interaction does not couple to the charge density since 
i s trace-less. Thus, the Thomas-Fermi screening 

for the coefficient A is absent. The local phonon Green 



function is D(u)i) — 2uqUJu/{i^>i 



'-,), which is given 



by the Fourier transformation of (T t m ai (t)m /j (0)) (/i = 
X, Y, Z) . uq = y/H/2RI-F C (jj£> is the position uncertainty 
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of Fe-ions, wd is the phonon frequency, and lui = 2-nlT 
is the boson Matsubara frequency. Then, for both Fe*- 1 - 1 
and Fe^ 2 ) , the phonon-mcdiatcd interaction is given by 

^24,42 - ¥34,43 = - {2 Aa? / if D {uj t ) = - ff ( W ,), 
¥22,22 = ¥33,33 = — ¥22,33 = _ 9\Ul), (2) 

as shown in Fig. 1 (b). Note that Vii',mm' is symmetric 
with respect to I -O- I', m -O- m! , and (ZZ') O (mm 1 ). We 

obtain g(0) 0.4 eV if we put i? Fc -A s ~ 2.4 A, a « 0.77 
A (Shannon crystal radius of Fe 2+ ), and wd ~ 0.018 eV. 
We have neglected the e-ph coupling due to d-p hybridiza- 
tion [14] considering the modest d-p hybridization in iron 
pnictides [15]. Thus, we obtain the multiorbital HH 
model for iron pnictides by combining eq. (2) with the 
on-site Coulomb interaction; the intra- (inter-) orbital 
Coulomb U (U'), Hund coupling J, and pair- hopping J'. 
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FIG. 1: (Color online) (a) FSs in the unfolded Brillouin 
zone, (b) Phonon-mediated electron-electron interaction, (c) 
A bubble-type diagram that induces the critical orbital fluc- 
tuations between (2,4) orbitals. (d) A ladder-type diagram 
that is ignorable when wd <C Ef- 

Now, we study the rich electronic properties real- 
ized in the multiorbital HH model [16]. The irre- 
ducible susceptibility in the five-orbital model is given 
by X?y, mm >(q) = -(T/N) E k GUk + q)G° mn ,{k), where 
G°(k) — [ie n + /i — -fffc] -1 is the ci-electron Green func- 
tion in the orbital basis: q — (q,ui), k — (k,e n ), and 
e n = (2n + 1)ttT is the fermion Matsubara frequency, fi 
is the chemical potential, and H% is the kinetic term given 
in Ref. [1]. Then, the susceptibilities for spin and charge 
sectors in the random-phase-approximation (RPA) are 
given as [17] 



(3) 



For the spin channel, j t t — U , U' , J, and J' for 
h = h = fa = fa, fa = fa ^ fa = fa, fa = fa ^ fa = '4, 
and Zi = I4 7^ fa = Z3, respectively [1]. For the charge 
channel, r c = —C — 2¥(a>;), where V{u>i) is given in eq. 
(2), and C h i 2< i aU = U, -U' + 2 J, 2ZJ' - J, and J' for 
Zi = fa = Z3 = Z4, Zi = Z 3 7^ fa — fa , fa = fa 7^ fa = fa, 
and h — h ^ h = fa, respectively [1]. Figure 1 (c) 



shows one of bubble diagrams for (2,4)-channel due to 
the "negative exchange coupling ¥24,42" that leads to a 
critical enhancement of x c (q) [18]. We neglect the ladder 
diagrams given by ¥(w;) in Fig. 1 (d) since wd -C Wband 
[8, 10]. We put lu d = 0.02 eV, U'/U = 0.69, J/U = 0.16 
and J = J', and fix the electron number n = 6.1 (10% 
electron doping); the density of states per spin is N(0) = 
0.66 [eV -1 ]. Numerical results are not sensitive to these 
parameters. We use 128 2 fc-meshes, and 512 Matsubara 
frequencies. Hereafter, the unit of energy is eV. 
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FIG. 2: (Color online) (a) Obtained U-g(0) phase diagram, 
(b) Obtained X24,42(<?,0) and X22,22(<7,0) for etc = 0.97. 

Figure 2 (a) shows the obtained U-g(0) phase diagram. 
a s (c) is the spin (charge) Stoner factor, given by the maxi- 
mum eigenvalue of T S ^X°{ ( 1, 0)- Then, the enhancement 
factor for x s ( c ) is (1 — a^^ 1 , and a s / c -\ — 1 gives the 
spin (orbital) order boundary. Due to the nesting of the 
FSs, the AFM fluctuation with Q ss (it, 0) develops as 
U increases, and s±-wave state is realized for a s < 1 [1]. 
In contrast, we find that the orbital fluctuations develop 
as g{0) increases. For U = 1, the critical value g C r(0) for 
a c = 1 is 0.4, and the critical e-ph coupling constant is 
A cr = g cr {0)N{0) = 0.26 [19]. Since the obtained A cr is 
close to A given by the first principle study [8] , strong or- 
bital fluctuations are expected to occur in iron pnictides. 
At fixed U, A cr decreases as J/U approaces zero. 

Figure 2 (b) and (c) show the obtained Xw mm 1 (9> 0) 
for {IV, mm!) = (24,42) and (22,22), respectively, for 
U = 1.14 and a c = 0.97 {g{0) = 0.40): Both of them 
are the most divergent channels for electron-doped cases. 
The enhancement of (24, 42)-channel is induced by the 
multiple scattering by ¥24,42- The largest broad peak 
around q = (0,0) originates from the forward scat- 
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tering in the electron-pocket (FS3 or 4) composed of 

2 ~ 4 orbitals. (FS1,2 are composed of only 2 and 

3 orbitals.) These ferro-orbital fluctuations would in- 
duce the softening of shear modulus [9], and also rein- 
force the ferro-orbital-ordered state below Ts [20] that 
had been explained by different theoretical approaches 
[21]: The divergence of X24,42 (X34.43) pushes the 2,4 
(3,4) orbitals away from the Fermi level, and the Fermi 
surfaces in the ordered state will be formed only by 3 
(2) orbital, consistently with ref. [20]. The lower peak 
around Q = (n, 0) comes from the nesting between hole- 
and electron-pockets. Also, the enhancement of (22, 22)- 
channel for Q = (tt, 0) is induced by the nesting via multi- 
ple scattering by V22,22 and V22,33- In contrast, the charge 
susceptibility J2i m Xu,mm( ( l^) 1S finite even if a c 1 

Since X22.33 ~ ~X22,22' 

Now, we will show that large orbital fluctuations, 
which are not considered in the first principle study of 
T c [8], can induce the s++-wave state when g(0) > 0. 
We analyze the following linearized Eliashberg equation 
using the RPA [1], by taking both the spin and orbital 
fluctuations into account on the same footing: 



k' ,m,i 

^^mim2(^ )^m2ni3 77147713 ( & )i(^0 



where W(q) = -^T s X s {q)T s + ±T c x c (q)T c - ±(r s - r c ) 
for singlet states. The eigenvalue Ae increases as T — s- 0, 
and it reaches unity at T = T c . In addition, we take 
the impurity effect into consideration since many iron 
pnictides show relatively large residual resistivity. Here, 
we assume the Fe site substitution, where the impu- 
rity potential I is diagonal in the d-orbital basis [3]. 
Then, the T-matrix in the normal state is given by 
T(e„) = [I- 1 - A r_1 EfcG(fc, en)]" 1 in the orbital basis 



[3]. Then, the normal self-energy is E n (e n ) 



T(e n ), 



where 



''imp 



is the impurity concentration. Also, the lin- 



earized anomalous self-energy is given by 

^ti'( e n) — j^ T P ^""^ Tj mi (e n )G mim , 2 (fc, e n ) A m2m3 (fc, e n ) 



X m/ x ra;i ( fc; ^■n) r ^l'm4,\ ^n)* 



(5) 



Then, the Eliashberg equation for 7ii mp ^ is given by 
using the full Green function G(k) = [ie n + fx — — 
£™(e„)] _1 in eqs. (4) and (5), and adding E^,(e n ) to 
the right hand side of eq. (4). Hereafter, we solve the 
equation at relatively high temperature T = 0.02 since 
the number of fc-meshes (128 2 ) is not enough for T < 
0.02, due to the fact that k-p in iron pnictides is only 1/5 
of that in cuprate superconductors. 

Figure 3 shows the ni mp -dependence of Ae at a c — 
0.98, for U = 1.11, 1.14 and 1.18. Considering large A E > 
0.8 at T = 0.02, relatively high-T c (< 0.02) is expected. 
For the smallest U (U = 1.11; a s — 0.85), we find that 
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FIG. 3: (Color online) nimp-dependence of Ae at a c = 0.98. 
If we put g(0) = (s±-state), Ae at rii mp = decreases by 
0.1 ~ 0.15, since the ferro-orbital fluctuations enhance both 
s++ and s± wave states. Inset: a c -dependence of Ae- 



nearly isotropic s++-wave state is realized; the obtained 
Ae is almost independent of ni mp , indicating the absence 
of impurity effect on the s++-wave state, as discussed in 
Refs. [3, 22]. For the largest U (U = 1.18; a s = 0.91), s±- 
wave state is realized at ni mp = 0; Ae decreases slowly as 
?iim P increases from zero, whereas it saturates for rij mp > 
0.05, indicating the smooth crossover from s±- to s++- 
wave states due to the interband impurity scattering. For 
U = 1.14 (a s — 0.88), the SC gap at n imp = is a hybrid 
of s_|__|_ and s±; only Aes2 is different in sign. 

The inset of Fig. 3 shows Ae for ,s++-wave state in 
the presence of impurities (nj mp = 0.15): Since Ae(o:c = 
0.98) — Ae(o; c = 0.90) is only ~ 0.15 for each value of 
U, we expect that relatively large T c for s ++ -wave state 
is realized even if orbital fluctuations are moderate. We 
stress that the obtained Ae is almost constant for w D = 
0.02 ~ 0.1, suggesting the absence of isotope effect in 
the s++-wave state due to the strong retardation effect 
[14]. By the same reason, Ae for the the s ++ -wave state is 
seldom changed if we put U = 3 in the Hartree-Fock term 
i(T s — r c ) in W(q), indicating that the Morel- Anderson 
pseudo-potential almost saturates. 

Here, we discuss the case U — 1.18 in detail: Figure 4 
shows the SC gap on the FSs in the band-representation 



for (a) n- u 



0, (b) 0.03, and (c) 0.08. They satisfy 



the condition N x J2k im |A; m (fc)| 2 = 1. The horizontal 
axis is the azimuth angle for the fc-point with the origin 
at T (M) point for FS1,2 (FS4); = corresponds to 
the fc^-direction. In case (a), s±-state with strong imbal- 
ance, [Apsil) |Afs2| <C Afs4, is realized, and Afs4 takes 
the largest value at 9 — tt/2, where the FS is mainly 
composed of orbital 4. In case (c), impurity-induced 
isotropic s++-statc [23] with Afsi ~ Afs2 ~ Afs34 is 
realized, consistently with many ARPES measurements 
[24]. In case (b), Aj- on FS1 is almost gapless. However, 
considering the fc z -dependence of the FSs, a (horizontal- 
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type) nodal structure is expected to appear on FSl, 2. 
In real compounds with T c ~ 50K, the s± — ¥ s++ 
crossover should be induced by small residual resistivity 
Pi m p ~ 20 /if2cm (ni mp ~ 0.01 for / = 1), as estimated in 
Ref. [3]. 
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FIG. 4: (Color online) SC gap functions for U = 1.18 as func- 
tions of 8 at (a) riimp = 0, (b) 0.03, and (c) 0.08, respectively. 



We comment that at n lmp = 0, s±-wave state is re- 
alized in the RPA even if a s < a c , due to factor 3 in 
front of \Y s x s (q)T s in W(q). For the same reason, how- 
ever, reduction in a s (or increment of U cl - for a s = 1) due 
to the "self-energy correction by {/" is larger, which will 
be unfavorable for the s±-wave state. Therefore, self- 
consistent calculation for the self-energy is required to 
discuss the value of a CiS and the true pairing state. 

Here, we discuss where in the a s -a c phase diagram in 
Fig. 2 (a) real compounds are located. Considering the 
weak T-dependence of \jT\T in electron-doped SC com- 
pounds [25], we expect that they belong to the area a c ^> 
a s . Then, s ++ -wave SC state will be realized without (or 
very low density) impurities, like the case of U = 1.11 
or 1.14 in Fig. 3. On the other hand, impurity-induced 
s± — > s ++ crossover may be realized in BaFe2(Asi_ x P x )2 
(undoped) or (Bai_ x Ka;)Fe2As2 (hole-doped) SC com- 
pounds, where AFM fluctuations are rather strong. 

Finally, we discuss the non-Fermi-liquid-like transport 
phenomena in iron pnictides. For example, the resis- 
tivity is nearly linear-in-T, and the Hall coefficient Rn 
increases at lower temperatures [4, 26]. Although the 
forward scattering induced by ferro-orbital fluctuations 
might be irrelevant, antiferro-orbital and AFM fluctua- 
tions with Q = (tt, 0) are expected to cause the anoma- 
lous transport, due to the current vertex correction [27]. 

In summary, we have proposed a mechanism of s++- 
wave SC state induced by orbital fluctuations, due to the 
phonon-mediated electron-electron interaction. Three 
orbitals (XZ, YZ and X 2 ~ Y 2 ) are necessary to leads the 
ferro-orbital fluctuations. The SC gap structure drasti- 
cally changes depending on parameters a s , a c , and n; mp , 
consistently with observed rich variety of the gap struc- 
ture that is a salient feature of iron pnictides. Orbital 
fluctuation mediated s ++ -wave state is also obtained for 



hole-doped cases, although the antiferro-orbital fluctua- 
tions becomes stronger than the ferro-orbital ones. 

The s-wave superconductivity induced by orbital fluc- 
tuations had been discussed in Ref. [17] for U' > U; this 
condition can be realized by including the Ai 9 -phonon 
[28]. In the present model, however, ^4i ff -phonon is neg- 
ligible since <? C r(0) given by ^4i s -phonon is much greater 
than <7cr(0) ~ 0.4 in Fig. 2 (a): The ferro-obtital fluctu- 
ations in Fig. 2 (b) originate from the negative exchange 
interaction caused by i? g -phonon, as shown in Fig. 1 (c). 

We thank D.S. Hirashima, M. Sato, Y. Matsuda, Y. 
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Note added in proof: After the acceptance of this work, 
we found that <? C r(0) ~ 0.4 in Fig. 2 (a) is reduced to half 
if all the e-ph matrix elements including 1, 5 orbitals are 
taken into account. Results similar to Fig. 3 are obtained 
by using g(0) ~ 0.2, whereas (vertical-type) nodes appear 
on FS3,4 during the s ++ — > s± crossover for U = 1.18. 
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